numerical integration with trapezoidal and simpsons rule